General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



FACILITY FOAM MI 


PREPRINT 


NASA TM X- 6S&36 

ACCELERATION OF PROTONS 
BY INTERPLANETARY SHOCKS 


KEITH W. OGILVIE 
JOHN F. ARENS 


FEBRUARY 1970 



GODDARD SPACE 

GREENBELT, 


N70-22675 


(NASA CH OA TMX OA AD NUMDERI 


• WUkl _ 

M 




/ 


K. W. Ogilvie 

Laboratory for Extraterrestrial Physics 
Goddard Space Flight Center 
Greenbelt, Maryland 


J. F. Arens 

Laboratory for Space Physics 
Goddard Space Flight Center 
Greenbelt, Maryland 


March 1970 


Space Plasma Physics Branch Preprint Series 


1 


Abstract 

Observations of charged particles at the time of passage of 
interplanetary shocks past the satellite Explorer 34 are discussed. 
The short duration increases in flux seen at 1 Mev are interpreted as 
particle acceleration, and are found to be consistent in duration 
and magnitude with the idea of energy gain by successive reflection 
between the earth's bow shock and the incoming propagating shock. 

The corresponding correlation length of the interplanetary field is 

_3 

deduced to be ~5xl0 A.U. From the small sample observed it appear^ 
that shocks occurring when the interplanetary magnetic field at the 
observers position does not intersect the bow shock do not show 
particle flux increases, and that particles from the solar wind are 
not accelerated to 1 Mev by any of the observed shocks. 


2 


Introduction 

Observations by Asbridge (1968), and by Frank (1967) have demon- 
strated the existence of particles streaming from the earth's bow 
shock in the general direction of the sun. These particles have 
energies greater by a factor of three or four than those characteristic 
of the solar wind protons. This might be interpreted as acceleration 
of particles of energy of order 1 KeV at the standing bow shock, and 
suggests the possibility of particle acceleration by interaction 
with propagating interplanetary shocks. 

Van Allen and Ness (1967) showed simultaneous observations of the 
magnetic field and of approx. 0.5 Mev protons on Explorer 33. At 
the time of passage of an interplanetary shock past the satellite, 
which was situated about 4xlo\m from the earth, a discontinuous drop 
in particle intensity was observed. As the particles originated at 
the sun simultaneously with the shock, the average particle on the 
outward side of the shock before it passed the spacecraft would be 
expected to have higher energy than the average particle on the sunward 
side, as a result of collisions with the shock. Such acceleration 
has been suggested by Vernov et al. (1969), in connection with so- 
called "Storm Particle Events". The increases discussed here, of 
duration about 1000 seconds, are not to be identified with such events. 
Axford and Reid (1963) have examined increases of low energy (~1 Mev) 
particle flux before the occurence of an ssc, interpreting these 
enhancements in terms of acceleration by successive reflections between 
the incoming shock and the earth's bow shock. 
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An orbit-theoretical treatment of the process of reflection of 
solar wind particles at the earth's bow shock has been published by 
Sonnerup (1969). Subject to the assumption that the particle energy 
is conserved in the coordinate frame in which the interplanetary electric 
field vanishes, acceleration of particles by about four times is 
predicted. The energising mechanism is the interplanetary electric 
field, doing work on the particles, supposedly trapped in the shock 
during the reflection process. This mechanism is discussed further below. 

Events showing apparent particle acceleration have also been 
studied by Singer (1970), who lists most of the events studied here 
and a number of others. He favors the hypothesis of the acceleration 
of particles in the shock front itself. A treatment cf the reflection 
of relativistic charged particles by plasma shocks, for the case 
where the shock thickness is less than the gyro-radius has been given 
by Hudson (1965). This condition applies to the present case. The simple 
non-relativistic acceleration mechanism discussed below assumes a high 
reflection coefficient. The magnetic field is not deviated though a 
large angle at the passage of an interplanetary shock, Ogilvie and 
Burlaga (1969). Under this condition, Hudson finds the reflection 
coefficient to approach unity for particles with large pitch angles, 

>75°. Unfortunately, pitch angle information was not available for the 
events discussed here. 

In this paper we combine observations of protons in the range 
1-10 Mev with plasma and magnetic field observations, and with the 
help of information on the solar proton spectra obtained on the same 
satellite and kindly supplied by Dr. Lanzerotti, we show: 
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1) . Shocks in the interplanetary medium do not always accelerate 
protons with solar wind energies to energies of order 1 Mev. 

2) . If particles with energies of a few hundred kev exist in the 
interplanetary medium at the time of arrival of a shock, an increase 
in flux above a threshold, (say 1 Mev) is sometimes observed. 

Increases produced by the mechanism discussed below can only be recorded 
when the particle spectrum is sufficiently steep and the interplanetaty 
magnetic field is suitably oriented. The flux drops to the ambient value or 
below after the passage of the shock. A mechanism for this effect, 

being essentially that proposed by Axford and Reid, will be discussed 
below. 

Experimental 

The identification of shocks and the measurement of their properties 
was carried out using the GSFC-University of Maryland plasma experiment 
and the GSFC magnetic field experiment on Explorer 34. Both of these 
instruments have been described before (Ogilvie and Burlaga, 1969) in 
a reference where the properties of most of the shocks used in this 
study are also described. The method of determination of the fluid 
quantities from the raw data is treated in Burlaga and Ogilvie (1968). 

The magnetic field observations were used for two purposes; to 
unequivocally identity whe passage of an interplanetary shock past 
the spacecraft,, and to establish the direction of the magnetic field 
vector at the shock front. The shock normal direction could sometimes he 

obtained, using the methods described in Ogilvie and Burlaga (1969). 



The l- 10 Mev protons were detected by an 86p. thick surface 
barrier detector which looked perpendicular to the satellite spin 
axis, the axis being normal to the ecliptic plane. The opening 
aperture was round with 60° full diameter. Each data sample 
accumulated counts from several revolutions, so azimuthal anisotropy 
information was not obtained. 

Results 

Between May 30, 1967 and 11 Jan. 1968 fourteen interplanetary 
shocks were detected by Explorer 34. Some details of these events are 
set out in Table I. Seven of them have been used as examples of 
hydiomagnetic shocks in the solar wind to test the applicability of 
the Rankine-Hugoniot relations; the remaining seven were not used, for 
a variety of reasons which are given in Table I. The observations made 
by the 1-10 Mev channel of the energetic particle apparatus were 
examined around the time of occurrence of all 14 events to look for 
evidence of particle acceleration or deceleration. 

In some events an increase in particle flux was observed starting 
about 20 minutes before the passage of the shock, and reaching a peak at 
the shock time, afterwards decaying to the pre-shock value 
or below. An increase in particle flux (nv) can be due to an increase 
in density or speed or both. Thus the mechanism must be identified 
in order to prove the occurrence of acceleration. Five events noted 
in Table I show appreciable flux increases. Two of these, which occurred 
on 30 May and 29 November 1967 are shown in Figures 1 and 2. The inset 
in Figure 1 shows the event on a larger time scale, and is included to 
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emphisi?;e the short duration of the increase. The quantity "Ambient 
rate" in Table I represents the counting rate which was observed in 
the 1-10 Mev channel immediately before the shock. It will be seen 
at once that all of the events showing acceleration occurred when 
this rate was &10, and that none of the shocks which occurred when 
the ambient rate was <10 showed an increase in the number of particles 
greater than 1 Mev at the time of the shock. It then appears that for 
these events the action of interplanetary shocks will not accelerate 
protons from the solar wind energies, (of order 1 KeV) , to energies 
of order 1 Mev contrary to the suggestion of Vernov (1969). In column 
eight of Table I we note whether the detection of the shock was 
immediately followed by a decrease in particle flux below the previous 
ambient level, of the type noted by Van Allen and Ness. A total of 


five such occurrences was seen. 
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The decrease referred to m column 3 is similar to that discussed by Van Allen and Ness. 
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The events on June 26, Aug. 11 and Sept 19 took place at times 
when the ambient flux was ^10, and yet rio appreciable flux increase 
at energies >1 Mev was observed. Events for which an increase was 
observed in the presence of an ambient solar proton flux are compared 
with these three in Table II. It will be seen that the only striking 
difference between the two classes of events in this small sample are 
the values of the ratio R. This quantity, the ratio of the fluxes in 
2 channels of the Bell Telephone Laboratory energetic particle experi- 
ment, is a measure of the steepness of the ambient solar proton spectrum. 
Its average value for the five events showing increases is 34 8, while 
that for the three events not showing increases is 9.6. 

Thus a steep spectrum is required for a large increase to be seen. 
This conclusion is exhibited . aantitatively in Fig. 3, where we see 
a plot of the expected increase in particle flux in a differential 
energy window against the differential spectrum exponent y. The 
corresponding values of the ratio \ are also shown, and points for 
the eight events in Table II are plotted. The diagonal lines represent 
the increase which would be observed as a furction of y, if the 
particles in this region of the spectrum had their energies increased 
by a.“l, 2, 3, etc. times. 

The events showing increases require, on this acceleration 
hypothesis , values of a of three or four. The three anomalous events, 

26 June, 11 August, and 19 Sept., are characterised by relatively flat 
spectra. Nevertheless, if the value of a had been 3 or 4 for these 
events. Increases of about a factor of two would have been seen, so 
an explanation for this effect must be sought . 
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Table II 



I max^ 1 before 

R 

0 

0 B 

H 

May 30 

20 

16 




June 5 

10 

41 

10° 

o 

-60 

316° 

Sept 20 

4 

10 

57° 

12° 

264° 

Nov 29 

5 

31 

47° 

10° 

151° 

Jan 11 

20 

76 

75° 

0.5° 

258° 

Events with ambient solar protons not showing 

acceleration. 

June 26 

1 

9 

30° 

-27° 

o 

r>* 

CM 

Aug 11 

1 

10 

31° 

21° 

300° 

Sept 19 

1 

10 

o 

00 

-24° 

272° 


Flux (0.56 - 0.60 Mev) 

~ Flux (1.2 - 2.4 Mev) 

®B» hourl y ave *age values. 


Note: - 
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III. Acceleration of particles by reflection between a moving shock 
and the bow shock. 

We assume the guiding centers of the particles (protons) to be 
constrained to move with speed V ^ along the interplanetary field lines. 

If the magnetic field line through the point of observation does not 
intersect both the bow shock and the moving shock we assume that 
there wi 7 l be no acceleration. The flux ratio is also a function of 
Y, the exponent of the differential energy spectrum of the particles, 
in agreement with the experimental results discussed above. 

We now estimate the increase in flux which would be seen by a 
detector situated at the point A in Figure 4, as the incoming shock 
mover with speed U along its normal direction. We assume that the particles 
pick up energy AE upon reflection at the moving shock and are mert- \y 
turned around at the bow shock. Thus the interval of time during 
which acceleration takes place is the time required for the point 
of intersection of the field line through the observers position with 
the shock to move from a point distant from the intersection with 
the bow shock to A, distant from the bow shock. The length L is 
to be identified with the correlation length characteristic of the 
motion of protons of this energy in the medium. A test of the 
hypothesis will be to determine that the necessary flux increase can 
be obtained with a reasonable value of L^. 

Let n be the normal to the moving shock and the angle between this 
vector and B be 0. The assumption that B intersects the bow shock 
is implicit. We assume that V x , the velocity of the particle 
perpendicular to the field line, is unchanged by the collision, but 
V ^ is increased by 2U Sec 0 by interaction with the shock. 
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1 , 2 2 2 2 
Thus, E , the energy after reflection, =%m(V j +V ^+4UV^SecB+4U 'Sec g) 

and this result does not hold for large values of g , 

E 1 -E , 4UV Secg + 4U 2 Sec 2 g 

i = U * 

Ei V 2 


4 U Secg»V,. , . „ -vvj.-v 

_ _JL (since V »U) 

V V 


= f.4 U Secg 

V 

where f is the cosine of the pitch angle. The sense of the field is 
unimportant 

Since the number of collisions per second with either shock is 
V/2L’ w **ere L *- s t ^ ie distance between the 2 shocks along the extended 
magnetic field line, 


AE = 2U -f - Sec g-At 
E L 


while the shock moves a distance UAt along its normal, its intersection 

with B moves a distance AL, so that UAt =-AL Cos g, and 

AE = -2 • f AL 
E L 

2f 

° r ^ =( H) 

Ei L 2 


Thus the effect of the motion of the shock from a remote point 


to the observers position is to increase the energy of each particle 

2f 

by a factor k= (Li / ) , which does not depend upon E . 

x, L2 1 

Suppose the particle spectrum has the form dN = AE ^ , so that 

dE 


the number of particles in a differential energy interval AE is 

dN AE, and N is the total number of particles. If all the particles 
dE 


have their energy multiplied by a constant factor k, then those which 

occupy AE at E after the acceleration are those which occupied AE at 

E^ before. Thus the ratio of the fluxes observed is 

fJ. 1 = V(E) A k^E^AE = k Y_1 , where A*(y-1)N. 
/F 1 V(E) A E Jy AE 
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This is the increase in flux due to acceleration; in addition there is 
an increase in differential density due to the compression alone, so 


that finally 


F 2/Fi = K 


Y-l 


= 


2 f(Y-l)+l 


( 2 ) 


Note that the increase according to this model would be the same 
at any point along L. As L decreases, however, the observed flux 
increases and should be a maximum at the time of shock passage, in 
agreement with observation. 

In Figure 5 we see a plot of the ecliptic plane with the average 
positions of the earth's bow shock and magnetopause superimposed. 

The position of the spacecraft and the directions of the magnetic 
field and shock normal on the ecliptic plane are shown for eight events. 
These are ones for which increases were seen (5) and also events 
for which an increase was not seen (3) despite the presence of an 
appreciable flux of solar protons before the event, Table II. The 
magnetic field vectors plotted are averages for the hour preceding 
the event, and the normals are those determined by the methods 
discussed in Ogilvie and Burlaga, 1969. For the events showing increases 
the magnetic field vector through the point of observation intersects 
the bow shock in each case for which we have data. For the 'anomalous' 
cases. Sept. 19 and Aug. 11 are clearly ones where B did not intersect 
the bow shock. Note that the value of the angle 0g does not affect 
this result. For June 26, the magnetic field was so disturbed that 
particles may not have been trapped on it. Thus we see that geometrically 
the idea of reflection between the two shocks is possible. 

We now determine whether the energy gains required are compatible 
with reasonable values of f and L^. For this purpose we can use the 
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September 20 and November 29 events, for which we know all the relevant 
quantities, except f. If Hudson's treatment is correct f is of order 0.1. 

TABLE III 


Event 

X ¥ 

h 2 

L 1 (f=l) 

L^f-0. 1) 

Sept 20 

3.3 2.03 

3.8xl0^cm 

~4x10" 3 A.U. 

-3 

~Url0 A.U. 

Nov 29 

4.5 3.1 

1. 3xl0^cra 

~2x10‘ 3 A.U, 

~4x10 _3 A.U. 


We identify this length with the correlation length of the interplanetary 
field. Thus when the moving shock is too far away from the bow shock, 
the particles are lost from the region before they make enough 
collisions to be accelerated. Since the values obtained are entirely 
reasonable it appears that all the features of this set of observations 
may be explained by multiple reflections between the bow shock and 
the advancing shock. 

We must now discuss the relationship of the mechanism proposed 
by Sonnerup, in which the energy of a solar wind particle is 
increased by four to six times at a single reflection at the bow 
shock, to the mechanism discussed here, where the energy gain per 
reflection is of order 10%. In the Sonnerup mechanism, a particle 
is quasi-trapped in the bow shock, where it can be accelerated by 
the action of the interplanetary electric field, which acts in a 
direction perpendicular to that of the shock motion. The order of 
magnitude of this effect for a particle of velocity V can be estimated 
by calculating the energy change due to accelerating a particle for 
a distance of one gyrodiameter in the interplanetary electric field. 



14 


E 2/ p _ 1 + 2eV sw B Sin ^ ( 2mv / eB ) £q 11 in Sonnerup (1969) 

b l 2 

rav 

~ * + ^ B in \Jr is the angle between the 

V 

magnetic field and the radial direction. 

9 -1 7 _i 

For a 1 MeV particle V = 1.4x10 cm sec , and V = 4x10 cm sec , 

sw 

Er ; ~ 1.1, but for a particle trapped in the shock so that V ~ V 
'E i “ sw 

the energy gain is a constant. If the particle does not become 
trapped in the bow shock it will gain no energy from the electric 
field, and merely be turned around by interacting with the materia? 
behind the shock. It will then gain energy during reflection by 
conservation of momentum, as assumed above. 

Duration of the increased particle flux . 

By measuring the duration of the particle increases we can 
estimate the shock speeds by using the relation. 

At = (Lj-l^) Cos 3 
U 

~ L, Cos 3 

■_i (3) 

u 

Both the Sept. 20 and Nov. 29 increases lasted 20 minutes +2 minutes. 

The corresponding values of U from equation 3 are 455 Km sec ^ for 
Sept. 20 and 170 Km sec for Nov. 29. A direct determination of U 
for Sept. 20 is not available, due to the uncertainty in post-shock 
plasma density. The directly determined value of Nov. 29 is 295 km sec 
(Ogilvie and Burlaga, 1969, also see Erratum, ibid, 1970). 
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These shock speeds are of the correct magnitude, and the disagreement 
on Nov. 29 Is within the rather large probable errors. 

We have shown that particles of solar wind energies are not 
normally accelerated by the action of interplanetary shocks to energies 
of order 1 Mev. Flux increases with durations of about 20 minutes are, 
however, observed, and their characteristics are consistent with 
acceleration by reflection between the travelling shock and the earth's 
bow shock, as suggested by Axford and Reid. This hypothesis gives 
reasonable values for the correlation length characterising the inter- 
planetary magnetic field, and is consistent with the measured shock 
speeds. It implies that acceleration will not be observed at the passage 
of shocks past a spacecraft far from the earth. A second mechanism 
might also operate to produce acceleration, but in that case one would 
expect to have observed an increase at the time of at least one of 
the three 'anomalous' cases referred to above. However a more extensive 
study over a wide energy range would be required to rule out such a 
possibility. 
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Fig. 1 
Fig. 2 
Fig. 3 


Fig. 4 
Fig. 5 


FIGURE CAPTIONS 

The increase observed at Explorer 34 in the flux of 
1-10 Mev protons on 30 May 1967. 

The increase observed at Explorer 34 in the flux of 
1-10 Mev protons on 29 Nov. 1967. 

The expected flux Increase in a differential energy 
interval plotted against the spectral exponent y. The 
values of the Ratio R, and points corresponding to eight 
events are plotted. The diagonal lines represent the 
increases which would be observed if the particles had 
their energies increased a«2, 3 etc. times. 

The geometry , not necessarily in the ecliptic plane, of 
the reflection of particles between the shocks. 

A plot of the ecliptic plane, showing the position of 
Explorer 34 at each of the events discussed. Boxes are 
drawn about the three anomalous events. The average 
positions of magnetopause and bow shock are also shown. 
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